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ABSTRACT 

The near infrared extinction powerlaw index (/3) and its uncertainty is derived from three 
different techniques based on star counts, colour excess and a combination of them. We have 
applied these methods to 2MASS data to determine maps of f3 and near infrared extinction 
of the small cloud IC 1 396 W. The combination of star counts and colour excess results in the 
most reliable method to determine /3. It is found that the use of the correct /3-map to transform 
colour excess values into extinction is fundamental for column density profile analysis of 
clouds. We describe how artificial photometric data, based on the model of stellar population 
synthesis of the Galaxy (Robin et al. (120031) ). can be used to estimate uncertainties and derive 
systematic effects of the extinction methods presented here. We find that all colour excess 
based extinction determination methods are subject to small but systematic offsets, which do 
not affect the star counting technique. These offsets occur since stars seen through a cloud do 
not represent the same population as stars in an extinction free control field. 

Key words: ISM: clouds, ISM: dust, ISM: extinction, ISM: structure 



1 INTRODUCTION 



Dust plays a major role in the dynamics of star formation. The de- 
termination of its properties and distribution is hence a major chal- 
lenge towards a better understanding. There are two widely used 
methods (star counts and colour excess) to determine the distribu- 
tion and column density of dust in the interstellar medium. 

The original idea o f usin g star counts to me asure dust extinc- 
tion dates back to Wolf ll923l) . Later Bok ll956l) improved the ac- 
curacy by implementing star counts at the completeness limit of the 
images. The difference between the measured and in trinsic colour 
of a star (colour excess) was used by Lada et al. J 1994) to determine 
the ex tinction. This approach was optimised by Lombardi & Alves 
(2001) from multi-band photometry. Recently Lombardi 120051) in- 
vestigated how the two methods can be optimally combined. 

Star counts gauge extinction at a given wavelength. The colour 
excess requires the knowledge of the reddening law to obtain an 
extinction measurement. Further the reddening law is required to 
transform these values into the widely used optical extinction (Ay). 
This transformation law is not known a priori. There is reasoning 
that a universal e xtincti on law is valid for the diffuse interstellar 
medium (Draine 120031) ). In translucent and dense regions, how- 
ever, processes like coagulation of dust grains and mantle growth 
can lead to a ch ange i n the dust propert ies (W hittet et al. 120011) : 
Cam bresy e t al. 1200 lh : del Burgo et al. 120031) : del Burgo & Lau- 
reijs J2005)), which could result in a different extinction law (Lar- 
son & Whittet (2005)). Ultraviolet and cosmic ray radiation has 
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also some influence in t he den sity and properties of dust grains (e.g. 
del Burgo & Cambresy 120061) ). Froebrich et al. 120041) showed that 
the near infrared extinction powerlaw index (/?) within dark clouds 
can vary significantly from the value obtained for the diffuse inter- 
stellar medium. 

The purpose of the present paper is the following. We de- 
scribe an easy-to-use and fast method to determine /3 from multi- 
filter broad-band observations (Sect.|2j- This method is applied in 
Sect.|3|to an example cloud (IC 1396 W). We discuss the implica- 
tions of usin g the determined /3 instead of the standard value of 
1 .85 (Draine 12003I) ). for density profile analyses of dark clouds. In 
Sect.|4|we present a technique to test the accuracy of extinction de- 
termination methods by artificial photometric catalogues based on 
the model of stellar population synthesis of the Galaxy by Robin et 
al. 12003T) . We finally discuss systematic uncertainties introduced 
into extinction maps when using colour excess methods (Sect. l4.4l . 
Our conclusions are then put forward in Sect.|5| 



2 DETERMINATION OF THE NIR EXTINCTION 
POWER LAW INDEX 

There are several possible methods to determine the near infrared 
extinction power law index (J3). In this section we describe tech- 
niques based on star counts, colour excess and a combination of 
both. The description is done in a general way, in order to adapt the 
procedure to different datasets. We discuss in detail the uncertain- 
ties of each method. 
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2.1 General notes 

Each method discussed below will determine as a function of 
two variables (u, v). Those can be extinction and/or colour excess. 
Hence we can write j3 = f(u,v). Each variable consists of N dif- 
ferent values (pixels in the extinction and/or colour excess maps). 
Following Taylor, we describe the difference of the individually 
measured values from the mean value as: 



n n , -\ 90 , .80 

0,-0 = (m - u)^f- + ( Vi - v)^f 
8u 8v 

A possible estimator of the variance of 0, defined as 
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can thus be expressed by the following equation: 
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(it, v). They can be estimated as the square of the noise in the input 
extinction and/or colour excess maps. The estimator of the covari- 
ance of the two input variables 
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can be determined from the extinction and/or colour excess maps as 
follows: multiply both images and determine the mean pixel value 
of the product image in an extinction free region. Note that in an 
extinction free field the mean values u and v are zero. 

In the next three subsections we introduce three different 
methods to determine the function = /(it, v) and the resulting 
estimator of the variance (a a). 



FigureQ (top-left) shows, as an example, the relative uncer- 
tainty A0/0 = . /<t|//3 versus for typical star counts in J- and 
H-band 2MASS data. The three different lines in Fig.Q] (top-left) 
denote extinction measurements Aj of 4, 6 and 10a. For 2MASS 
data in the area of IC 1396 W we have determined a sc = 1.4 and 
7sc = 1.4 (see below). It is found that a high accuracy for the ex- 
tinction measurements is required to obtain a good signal to noise 
ratio (S/N) for 0; for typical values of = 1.4 at least a 4a detec- 
tion of the extinction in the J-band is needed to determine with 
an accuracy of 40 %. Note that this is valid for a single measure- 
ment only (1 independent pixel in the extinction map). The S/N can 
be improved by averaging over several independent measurements 
(see below). 

2.3 Colour excess 

The most simple approach of a colour excess method is to measure 
the apparent brightness of a star at two different wavelengths, ui\ 1 
and raj 2 , and compare the colour m\ 1 - m\ 2 with the theoretical 
value for main sequence stars or the colour of stars measur ed in a n 



extinction free control field mx\ - m" 2 (e.g. Lada et al. 119941) ) 
The colour excess (Ai — A2) is defined as: 



(Ai - A 2 ) = {m\ 1 - m\ 2 ) - (m/j - m/J 



(7) 



If we assume that this colour excess is entirely due to extinc- 
tion and, as for the star count analysis, a constant is valid, we can 
determine the extinction at A 2 by 



A x „ = 



(Ai - A 2 > 



(8) 
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Measuring (A2 — A3}, with Ai < A2 < A3, we as well can 
determine the extinction at A2: 



2.2 Star Counts 

Extinction mapping using accumulative star counts allows us to de- 
termine the extinc tion at wavelength A whe re the photometry is ob- 
tained (e.g. Wolf Jl923h : Cambresy et al. <2002l) : Froebrich et al. 
120041) ). It is assumed that extinction has been measured with two 
different filterbands at reference wavelengths (Ai < A2). Also we 
assume the extinction has a power-law dependence A\ oc A - ^ in 
the wavelength range of interest. Then, the unknown index can 
be determined from: 







\n{A x JAx 2 ) 
ln(A 2 /Ai) ' 



(5) 



To estimate the uncertainty of we introduce two parameters: 
1) a sc , which corresponds to the ratio of the estimators of the vari- 



ance in the two star count maps, i.e. q s . 
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which corresponds to the ratio of the estimator of the variance in 
the star count map at Ai and the estimator of the covariance of the 
two star count maps, i.e. 7 SC = ct\ x /(Ja x a x ■ The estimator of 
the variance of can then be expressed as: 
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The right hand side of both equations should essentially give 



the same value. We can define the colour excess ratio R : 
and obtain the following equation for 0: 
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Fig.|2| shows the solution of Eg . 11 01 for some filter combina- 
tions. These plots allow us to estimate for a certain filterband 
combination and measured colour excess ratio. 

Since cannot be determined analytically from Eg . 11 01 we 
assume that — B(i?) provides a solution of Eg. llOl with sufficient 
accuracy (the function B(i?) could be a polynomial of order n). 
Then we can determine the estimator of the variance of using: 
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Note that the partial derivative 8T5(R)/8R can directly be 

obtained from Eg . 11 01 by using the inverse function theorem. We 



define the parameters a co = ^(x 1 -x 2 )/ a2 (\2-^3) an ^ ^ c 
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Figure 1. Top-left: Relative uncertainty of /3 obtained from J- and H-band star counts. Top-right: Relative uncertainty of /3 obtained from (J — H) and 
(H — K) colour excess ratio. Bottom: Relative uncertainty of /3 obtained from the combination of J-band star counts and (J — H) colour excess. The plots 
show A/3//3 referring to measurements of Aj (left panels) and (J — H) (right panels) with a S/N of 4 (solid line), 6 (dotted line), and 10 (dashed line). Values 
of a BC = 1.4, 7sc = 1.4, ctcosc = 13, 7 CO sc = 6.5, a co = 3.6 and 7 co = 3.7 were used for the determination (see text for details). The two vertical lines 
in the plots indicate values of (3 = 1.4 and 1.8, which represent the value for dense regions and the diffuse ISM, respectively. Note that the two bottom panels 
are identical despite showing apparently different behaviour. The relative uncertainty of as a function of the S/N in the Aj map (Eg . 1131 can be converted 
into the relative uncertainty of ft as a function of the S/N in the (J — H) map using Eq J12l and the values of o C osc and 7 CO sc- 



<A 1 -A 2 )/ cr (A 1 -A 2 )(A 2 -A3) 



Average values based on 2MASS 
data for IC1396W are a co = 3.6 and j co = 3.7. Note that 
these values depend on the intrinsic scatter of the stellar colours, 
the spacial distribution of stars and the photometric uncertainties 
in the data, and are hence subject to change slightly with position. 
In Fig.Q (top-right) the resulting relative uncertainty of (5 for this 
example is plotted. Similar to the star counting technique, high S/N 
colour excess maps are required to reliably estimate j3. We note, 
however, that for the same dataset usually a higher S/N is achieved 
for the colour excess method than for the star counting (see e.g. 
Cambresy et al <2002l) ). 



2.4 Combining star counts and colour excess 

If observations at only two wavelengths Ai < A2 are available, 
we do not necessarily need to determine /3 solely from star counts. 
We can combine star counts and colour excess. If we measure the 
extinction A\ 1 by star counts, the colour excess (Ai — A2), and 



assume that j3 is constant in the wavelength range of interest, then 
(3 can be determined as: 



/3=-ln 1- 



(Ai - A 2 ) 
Ax, 



/MA2/A1 



(12) 



As for the other two methods, we can estimate the uncertainty 
associated to the procedure. Similar to those we define the param- 
eters ttcosc = ffi Al /^(Ai-A 2 > and 7cosc = <ri Ai /<7> Xi<Al _A 2 >, 
and obtain the estimator of the variance of j3: 
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FigureQ (bottom) shows two examples of the relative uncer- 
tainty of P when determined from the combination of star counts 
and colour excess. Both cases combine J-band star counts and 
(J — H) colour excess. In Fig.Q (bottom left) J-band star counts 
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K 




Figure 2. Solution of Eg J 101 for a number of common filter/wavelength 
combinations. This plot allows us to estimate /3 from the ratio of the colour 
excess of the two shorter and the two longer wavelength. We plot: BVR 
(solid), VRI (dotted), RIJ (short-dash), IJH (dash-dot) and JHK (long-dash). 



are used as a measure of the S/N. In Fig.Q (bottom right) the 
(J — H) colour excess is used. Values ofa cosc = 13 and 7 COS c = 
6.5 are obtained for the example (IC 1396 W). 



2.5 Other combinations 

If observations at three wavelengths are available and we assume 
that there is a constant j3, in principle a variety of other possible 
combinations of colour excess and star counts emerges. One can 
combine colour excess between two of the wavelengths with star 
counts at a third wavelength. Also two other combinations using 
purely colour excess are possible. Since the procedure for the de- 
termination of (3 and its uncertainty, as well as the results are very 
similar to the techniques described in the above paragraphs, we re- 
frain from presenting them here. 



2.6 Usage of broad-band filters 

The analysis so far is implicitly carried out using the simplifying 
assumption that extinction and colour excess measurements are car- 
ried out at monochromatic wavelengths. Usually the observations 
are, however, obtained using broadband filters. Does this influence 
our results? 

Firstly we note that extinction and colour excess measure- 
ments are obtained by averaging data over a number stars in a cer- 
tain area. For example the colour excess maps of IC 1396 W (pre- 
sented below) are determined using on average 25 stars at each po- 
sition. The average spectral energy distribution of these stars indeed 
will change the reference wavelength of the used filter. For NIR ob- 
servations, as used here, the reference wavelength will be shorter as 
the reference wavelength of the used filter. Second we note that all 
determination procedures for j3 and its uncertainty depend on the 
ratio of two wavelength. A change of this ratio would imply sys- 
tematic off-sets in our determined /3- value. Since for all NIR filters 
the shift in reference wavelength is toward shorter wavelength, the 
resulting change in the wavelength ratio will be small. 

To estimate this we have convolved black body curves with 
temperatures between 2000 K and 12000 K with 2MASS filter pro- 
files. It is found that the resulting shift in the reference wavelength 
ratio is at most 0.5 %. This generally results in about 3 % system- 
atic off-sets for the value of j3. It is hence much smaller than the 



determined statistical noise, which is in the order of 30 % or higher. 
Thus, the use of broad band filters does not influence our results. 



2.7 Discussion 

In general all methods described here have similar uncertainties 
when applied to the same wavelength range. But it is possible to 
determine the method leading to the highest S/N in /3 for any par- 
ticular dataset. However, there might be additional systematic un- 
certainties (see Sect. 14.41 for a discussion). FigureQ allows us to 
compare the different methods. In the left panels the relative un- 
certainties of j3 are plotted as a function of the S/N in the J-band 
star count map. Note that for purely star counting the uncertain- 
ties are larger as for the combination of star counts and colour ex- 
cess. In the right panels of Fig.Qthe relative uncertainties of /3 as a 
function of the S/N in the (J — H) colour excess map are plotted. 
We find that the scatter is smaller when combining colour excess 
and star counts instead of using the colour excess ratio. This is es- 
pecially the case towards smaller values of j3. We conclude that 
for our example (IC 1 396 W) the combination of J-band star counts 
and (J — H) colour excess provides the best /3-map. Note that this 
might change for other datasets, since the estimator of the variance 
of j3 depends on the values of Qi and 71. 

Even when using the best method to determine /3, the uncer- 
tainties based on individual measurements (one independent pixel 
in the extinction/colour excess maps) are rather large even for high 
S/N extinction data. However, several pixels may be averaged to 
decrease the uncertainties. Although a simple averaging over ad- 
joining pixels would decrease the spatial resolution it is in general 
a good approach in regions were no apparent changes of /? are seen 
in the data (i.e. the scatter of /3-values is consistent with the statis- 
tical expectations; 68 % of the data points are within the la error 
bars, etc.). If a dependency of /3 on the column density is found 
(as for IC 1396 W, see below), we suggest to apply the following 
method to determine how (3 depends on the column density within 
individual clouds. 

• Accumulated star counts or its combination with colour 

excess: For each independent pixel /3 is plotted against extinction 
(e.g. A\ x ). A running mean in this diagram can then be used as rep- 
resentation of the relation P-A\ 1 . Alternatively mean values of /3 in 
certain ^Ai-bins can be fit by an appropriate function to establish 
the relation /3-Ax 1 ■ 

• Colour excess: The determination of the extinction from 
colour excess requires the knowledge of /3. Thus, as a first step 
for each independent pixel the value of /3 is plotted against the 
colour excess (e.g. (Ai — A2)). Similar to the accumulated star 
count method we establish the relation /3-(Ai — A2) using a run- 
ning mean or by fitting an appropriate function. This relation then 
can be used to determine the extinction from the colour excess and 
to determine e.g. the relation P-A\ 1 . 



3 A CASE STUDY: IC 1396 W 

In order to illustrate the /^-determination procedure and its impli- 
cations on extinction and cloud column density profiles, we present 
in this section an example, the small cloud IC 1396 W. This object 
is part of the Cep OB 2 association, about 8' in s ize and at roughly 
750 pc distance (e.g. Froebrich & Scholz (2003)). The position of 
the cloud in the sky and the 2MASS detection limits facilitate our 
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Figure 3. Top: Extinction maps of IC 1 396 W obtained from accumulative star counts in 2MASS data. Contours in each panel start at 0. 1 5 mag of extinction in 
the J-band (left), H-band (middle) and K-band (right) and increase by 0.15 mag with each step. Middle: Extinction maps of IC 1 396 W obtained from (J — H) 
colour excess in 2MASS data. Contours start at 0.15 mag of extinction in the J-band and increase by 0.15 mag with each step. The conversion of colour excess 
into extinction uses /3 = 1.8 (left), the /3-map obtained from the (J — H) and (H — K) colour excess ratio (middle) and the best /3-map determined from 
(J — H) colour excess and J-band star counts. Bottom: The same as the middle row but for (H — K) colour excess. Contours are in H-band extinction. 



analysis in the way that there is only a negligible fraction of fore- 
ground stars to the cloud. Hence we do not need to remove fore- 
ground stars when determining extinction maps. 



3.1 Extinction maps 

In the top panel of Fig.|3|we show the three extinction maps {A}, 
Ah, Ak) for this cloud obtained from accumulative star counts. 
The contours start at 0.15 mag extinction at the filter wavelength 
and increase by this amount. The basic structure (central peak and 
small extension to the East) of the cloud is the same in all three 
filters. The extinction maps obtained from the colour excess maps 
show a slightly different structure (see bottom panel of Fig.|3j- In 
particular we can trace much better the outer areas of the cloud. 
There is an apparent small offset between the star count and colour 
excess extinction maps. This can be explained by possible differ- 
ences in the selection of the control field (for the colour excess an 
apparently IRAS emission free field south of the cloud was chosen, 
in contrast to the surrounding l°xl° field for the star counts). 



3.2 /3-maps 

In the top panel of Fig.|4]we show the /3-maps of the cloud, as 
determined from the star counts in JH, HK and JK and using Eq.[5] 
The bottom left panel of Fig.|4]displays the /3-map obtained from 
the (J — H) and (H — K) colour excess maps and Eq.^j| The 
two remaining maps are determined using the combination of star 
counts, colour excess and Eg . 11 21 In the bottom middle panel we 
combine J-band star counts with (J — H) colour excess, and in the 
bottom right panel H-band star counts and (H — K) colour excess. 
In all cases only pixels are shown, where the used extinction or 
colour excess is 3a above the noise. 

The first thing to note is that /3-maps obtained using colour 
excess cover a larger area, indicating the superiority of the colour 
excess method to the star count method in the NIR, in terms of S/N 
in the extinction maps. We also see a structure in the /3-map that 
agrees well with the structure in the extinction maps. This suggests 
a dependence of /? on the column density in this cloud. In the outer 
regions of the cloud values for j3 close t o the canonical value for 
the interstellar medium of 1.85 (Draine 120031) ") are found. In the 
centre of the cloud, where extinction is higher, j3 drops to signifi- 
cantly smaller values, partly even below 1 .0. When using solely star 
counts to determine j3 a smaller area is above the 3er noise level and 
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Figure 4. Top: Gray scale /3-maps obtained from accumulative star count extinction maps for the wavelength range JH (left), HK (middle) and JK (right). 
Bottom left: Gray scale /3-map obtained from (J — H) and (H — K) colour excess maps. Bottom middle: /3-map obtained from the combination of 
(J — H) colour excess and J-band star counts. Bottom right: Gray scale /3-map obtained from the combination of (H — K) colour excess and H-band star 
counts. In all panels only areas with extinction or colour excess values above the 3<r noise level are used. Isocontours of /3=1.0 and 1.5 are marked on gray 
scale maps where darker colours represent higher values of /3. 



no structural information can be obtained. However, the combina- 
tion of J-band star counts and (J — H) colour excess agrees very 
well with the /3-map obtained solely from colour excess. The small 
offset in the absolute /3-values is due the different control fields 
used (see discussion above). We note, however, that differences in 
/3 within a single cloud can be easily determined by our method. 
The correct absolute value for /3 requires a carefully chosen com- 
parison field, which is free of extinction. 



3.3 Column density profile 

How does the change of /3 influence column density profiles of the 
cloud? We plot in Fig.[3]fhe column density profile (normalised to 
the peak) from East to West across the cloud at a declination of 
57.93°. The solid line shows the profile from J-band star counts 
and the dotted line from (J — H) colour excess (converted with a 
constant /3-value). There is a significant difference between the two 
profiles. The small peak in the East shows a much higher column 
density relative to the main peak when using colour excess instead 
of star counts. Note that this difference is unchanged when varying 
/3, as long as a constant value is used for the whole cloud. 

We used our /3-map obtained from the combination of J-band 
star counts and (J — H) colour excess to convert the (J — H)- 
map into column density. The resulting profile is shown in Fig.|5| 
as dot-dashed line. It is obvious that this dramatically changes the 
column density profile which now much more closely matches the 
profile obtained solely from star counts. Note that we obtain the 
same result when using the /3-map obtained from colour excess 
only (dashed line). 

This example shows that the column density structure of 
clouds can differ significantly depending on the method used to 



determine the extinction maps. This in particularly holds when 
changes of /3 are found within the cloud. Thus, extinction maps 
from accumulative star counts or colour excess (using the cor- 
rect /3 values) should be favoured when investigating the struc- 
tures/profiles of the denser parts of clouds, where changes of /3 
are more frequently observed. When investigating the outer areas 
of clouds, where /3 is very close to the ISM value, or at least con- 
stant, colour excess extinction maps provide the better solution due 
to their higher signal to noise. 



4 TEST OF EXTINCTION DETERMINATION 

Given the possible large uncertainties in the determination of /3, it is 
worthy to know the uncertainties of the underlying extinction maps 
from star counts or colour excess determination procedures. In par- 
ticular any systematic effect should be investigated and corrected 
for. Here we introduce a scheme to test the reliability of the extinc- 
tion determination methods. A detailed discussion about the influ- 
ence of the box size to the measured extinction, systematic offsets 
in colour excess methods and limitations for reliable foreground 
star detection is given. 

4.1 Artificial data 

We used the mo del of stellar population synthesis of the Galaxy 1 
by Robin et al. y003j) to generate artificial data. This model al- 
lows us to generate artificial photometric catalogues for a given 
region in the sky. Constraints for the limiting magnitude in each 

1 available at http://bison.obs-besancon.fr/modele/ 
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Figure 5. Column density profiles (normalised to zero at the eastern edge 
of the cloud and unity at the peak) across IC 1 396 W at a declination of 
57.93°. The solid line shows the profile obtained from J-band star counts 
and the dotted line is based on (J — H) colour excess and a constant /3. 
The dash-dot line uses a combination of J-band star counts and (J — H) 
colour excess to determine the /3-map, required to convert (J — H) into 
extinction. The dashed line is the column density profile when using the 
/3-map obtained from the (J — H) and (H — K) colour excess ratio. 
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Figure 6. Construction of the test field. An artificial star catalogue without 
extinction is generated for the whole 2°x2° area. Within a 0.5°x0.5° field 
(solid lines) this catalogue is exchanged by an artificial catalogue of stars 
with a cloud at a distance d and an extinction Ay. The final extinction map 
is determined for the central l°xl° area (dashed lines). We measured the 
extinction in the central 0.25°x0.25° field (dotted lines). 



filterband can be applied. We performed all our simulations us- 
ing 2MASS catalogue magnitude limits (J max = 16.5, H max = 16.0, 
K max = 15.5). As position we selected (£ = 81°, 6 = 0°). This coin- 
cides with the position of the DR 2 1 star forming region and hence 
the quantitative results of this study can be applied for this area. For 
any other region in the sky only the qualitative part of the analysis 
can be applied; for quantitative results the simulations described 
here have to be applied to the specific area. 

As mentioned, we determine only the extinction at a certain 
wavelength when using accumulative star counts, or we can deter- 
mine the colour excess between two wavelengths. The conversion 
into optical extinction requires the knowledge of the extinction law. 
Since we kno w the extinction law applied to generate the artificial 
data (Mathis ll990T) V we are able to directly compute the optical 
extinction. This is what finally is used to estimate the reliability of 
our extinction determination procedure. 

In Figure|6] we show the layout of our test field, used in the 
simulations. The following procedure was applied to generate arti- 
ficial observations for this 2°x2° field: 

(i) A photometric catalogue (set of stars) for a one square de- 
gree field around the coordinates (/ = 81°, b = 0°) was downloaded 
using the 2MASS magnitude limits, a standard galactic extinction 
law of Ay =0.7magkpc _1 and no clouds. To fill the 2°x2° field 
with this set of stars, each star was supplied with four different, arti- 
ficial, uniformly distributed coordinates, that placed it in the 2°x2° 
field. In other words, these artificial coordinates are pairs of ran- 
dom numbers generated to mimic a uniform distribution of stars in 
the sky. 

(ii) A photometric catalogue with the same magnitude limits and 
position was downloaded, but this time including a cloud with a 
certain distance and a certain extinction. Each star in this second 
catalogue was supplied with artificial, uniformly distributed coor- 
dinates, that placed it in the central l°xl° field. 

(iii) The two catalogues where merged in a way that all stars in 
the first catalogue which fell in the central 0.5°x0.5° area (solid 
lines in Fig.[6} are replaced by the stars of the second catalogue in 



this area. In other words we created an artificial catalogue, simulat- 
ing a squared 0.5°x0.5° cloud in the centre of a 2°x2° field. 

For each of these artificial datasets we determined the extinc- 
tion map for the central l°xl° field and measured the mean and 
median extinction of the cloud in the central 0.25°x0.25° . This set- 
up allows us to perform different tests to determine the reliability 
of our extinction determination method. 



4.2 Zero extinction (i.e. no clouds) 

Here we consider an object with zero extinction at zero distance to 
estimate the contribution of the simulated photometric datasets and 
artificial coordinates to the noise of the method. Such an object is 
equivalent to consider "no cloud", but it is also representative of a 
region with negligible extinction. 

We first fix the set of stars in the control field, the cloud and 
the coordinates assigned to these stars. This allows us to determine 
how the chosen box size influences the noise in the extinction map. 
As expected we find that the box size is inversely proportional to 
the noise, with very small deviations. This can be explained since 
the noise should be inversely proportional to the square root of the 
number of stars in each box, and then inversely proportional to the 
box size. As box size for all subsequent analysis we chose 1.2', 
which ensures on average 25 stars per box. 

We have performed several tests in order to disentangle the 
different components contributing to the noise in the extinction 
maps. In particular: a) the artificial coordinates of the stars in the 
field; b) the set of stars in the control field; c) the set of stars in the 
area of the cloud. The contribution to the noise from each of those 
three components can be determined from the scatter in the syn- 
thetic extinction maps, by varying one component and fixing the 
other two. 

Our results are summarised in TableQ There we compare the 
noise (i n mag nitudes of Ay, converted using the extinction law of 
Mathis (1990J)) between K-band star counts and (H — K) colour 
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Table 1. Contribution of artificial coordinates and set of stars on the noise 
in the extinction maps from star counts and colour excess for model clouds 
with negligible extinction. We show here results for K-band star counts and 
(H — K) colour excess methods. The noise is given in magnitudes of Ay. 



influences 



noise from 
star counts colour excess 



artificial coordinates 


0.15 


0.03 


set of stars in 






control field 


0.13 


0.05 


cloud 


0.10 


0.04 



total noise 



0.22 



0.07 



excess. We find that the star count method leads to three times 
larger noise (for our chosen magnitude limits and box size). When 
using the colour excess method, the set of stars in the control field 
or the cloud is responsible for most of the scatter. For star counts, 
however, the artificial coordinates dominate the noise. 



4.3 Clouds at zero distance 

Since we have established the influence of the simulated star sets 
and artificial coordinates, we now investigate the changes for dif- 
ferent cloud extinction values (Ay = 0.0, 0.2, 0.5, 1.0, 3.0, 5.0, 
7.0, 10.0 mag). To avoid uncertainties due to foreground stars, the 
model clouds are placed at zero distance. 

We test how the scatter in the extinction map changes with 
Ay. FigureQshows how the noise in the extinction maps changes 
with Ay of the model cloud. The dashed line corresponds to the 
colour excess method and the dotted line to the accumulated star 
counts. For low extinction values, below one mag of Ay, the mea- 
sured values agree very well with the estimated values in the above 
section. This can be explained since such low extinction values only 
marginally influence the number of stars in each box and therefore 
the noise in the maps. For larger Ay-values, however, the scatter 
increases significantly. This is largely due to the fact that a much 
lower number of stars is present in each box. When comparing the 
two methods we find that accumulated star counts lead on average 
to a three to five times higher noise. 



4.4 Systematic offsets — distance correction 

We simulated clouds with different extinction values at zero dis- 
tance. This ensures that no foreground stars are present and we 
should be able to fully recover the applied extinction values of our 
artificial clouds, given the fact that we kn ow the reddening law ap- 
plied to the artificial data (Mathis (l99fj)). Using accumulated star 
counts the measured extinction is in very good agreement with the 
real (input) value. The measured values scatter only 0.4 a around 
the input value and hence no measurable systematic effects exist. 
In absolute terms these deviations are larger than what is obtained 
by the colour excess method, due to the in generally higher noise 
of the accumulated star count technique. The offsets obtained by 
the colour excess method are hence small but systematic, i.e. above 
the internal scatter of the method. Fig.[8|shows the variation of the 
offset with Ay. In the following we discuss their cause. 

For small extinction values (Ay < 1 mag) the deviations 
are smaller than 1 a. For higher extinction (Ay up to 7 mag) we 
systematically underestimate the extinction by up to 2a. For very 
high column densities (Ay ~ 10 mag) this effect turns and we 





1.00 


1 i 1 i 1 i 1 i 1 i 1 


- 






/ 








/ 




> 


0.75 


/ 


- 




/ 










- 


C5 




/ 




a 


0.50 


/ 


- 


OJ 




/ 




or; 




/ 


- 


O 




/ 




2 


0.25 


...«■>• _^^ 


- 












4 6 

A v [mag] 



10 



Figure 7. Noise in the extinction maps as a function of Av of our model 
clouds. The solid line results from the colour excess method (H — K) and 
the dashed line from accumulated star counts (K-band). The two dotted 
horizontal lines mark the noise due to the set of simulated stars and artificial 
coordinates determined by our test in Sect. l4.2l for the two methods (the 
upper line is for the K-band star counts, the lower for the ) (H — K) colour 
excess. 
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Figure 8. Systematic offsets (solid line) of the recovered extinction values 
for our model clouds depending on the optical extinction when using the 
colour excess method. The dashed lines mark the ±lcr uncertainties. 



overestimate the extinction. Clearly an explanation for this effect 
is needed. Due to the setup of our test, a contribution from fore- 
ground stars can be neglected. For small values of extinction (be- 
low 1 mag Ay) the effect is negligible. This implies the reddening 
of the stars due to the dust in the cloud itself is the cause of the off- 
sets. When comparing the mean/median colour of stars in the cloud 
and an extinction free control-field, we have to keep in mind that in 
both areas stars down to a certain limiting magnitude are detected. 
We see, however, in these two fields different populations of stars. 
Stars behind the cloud, detected down to the same limiting mag- 
nitude, are on average at a different distance and redder. Thus the 
mean/median intrinsic colour of stars behind the cloud differs from 
the mean/median intrinsic colour of stars in the control field. In 
conclusion all methods using colour excess to estimate extinction 
are subject to this effect. Even if small for low extinction values it 
might become significant and systematic for denser regions, chang- 
ing effectively the clouds column density profile. 

The particular quantitative behaviour of the offsets is a com- 



Extinction techniques and impact on dust property determination 9 



plex interplay of various parameters. The position in the sky and the 
limiting magnitude of the observations determine the population 
of stars contained in the control field. The extinction of the cloud 
then determines how this population is changed when seen through 
the cloud. Extinction could have the effect that less red-giant back- 
ground stars are seen through the cloud. Also it is possible that 
less faint stars close behind the cloud are visible. If one of those 
effects prevails, systematic offsets are expected. The extinction fur- 
ther could cause the median distance of the stars seen through the 
cloud be different from the the stars in the control field. This can 
also introduce systematic offsets due to interstellar extinction (typ- 
ically 0.7 maeAv/kpc, but reaching up to 2 mag/kpc in certain re- 
gions; Joshi J2005J)). All these effects are further influenced by the 
actual distance of the cloud. Due to this complexity it is not possible 
to describe the offsets using a simple model with a few parameters. 
The only way to determine the quantitative offsets and to correct 
them is to perform simulations similar to the ones presented here, 
using the limiting magnitude and coordinates of the real observa- 
tions as well as the distance and extinction range for the observed 
cloud. We note that the quantitative effect most likely depends on 
the particular model used. However, it is important to realise its ex- 
istence and to know its possible order of magnitude. Thus it should 
be considered in the analyses of high extinction regions with colour 
excess methods. 

Do those offsets have a systematic influence on the determi- 
nation of (31 One can estimate that the value of j3, determined with 
or without the correction of the systematic offset, differs by at most 
5 %. This is usually much smaller than the scatter of/? (in the order 
of 30 % or higher) due to the scatter in the colour excess maps (see 
e.g. Fig.0. Only for very high S/N data with a very small internal 
scatter of j3 these offsets become important. 

4.5 Clouds with variable distance 

Foreground stars have to be identified and removed from the input 
catalogue because they influence the correct extinction determina- 
tion. This is usually done by identifying all stars in the field of the 
cloud that are much bluer than what would be expected from the 
extinction of the cloud. The notion is stars behind the cloud are 
red, stars in front of the cloud are blue. This requires, however, the 
detection of a significant extinction towards the cloud prior to the 
foreground star identification. Here we investigate the difficulties 
for the identification of foreground stars for the star counting and 
colour excess method. 

We analyse how the extinction of our model cloud depends on 
the distance of the cloud. Distances of 0,0. 1,0.5, 1, 2, 3,4and5kpc 
were considered. In Fig.|5|we plot how the measured extinction de- 
pends on the cloud distance for a 10 and 3 mag Ay cloud. Accu- 
mulative star counts (solid lines) and colour excess (dashed lines) 
extinction maps, prior to foreground stars selection, are shown. 

There are obvious differences between the two methods. The 
star counting method shows a steady decline of the measured 
extinction with increasing cloud distance. For the colour excess 
method at first only a very small decline with cloud distance is 
observed. After a certain critical distance, however, the measured 
extinction drops to zero. This is the distance at which more than 
half of the stars in the field are foreground to the cloud, and the me- 
dian colour represent foreground stars. The higher the extinction in 
the cloud the smaller is the maximum distance up to which a cloud 
can be detected using colour excess since less background stars are 
detected in the field. 

When investigating distant clouds, it is hence of advantage to 
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Figure 9. Recovered extinction depending on the distance of the cloud us- 
ing accumulative star counts (solid) and colour excess (dashed) methods for 
clouds with an optical extinction of 10 (lines stalling at 10 mag for zero dis- 
tance) and 3 mag (lines starting at 3 mag for zero distance). The horizontal 
lines mark the input cloud extinction. 

use the star count extinction maps for foreground star identification. 
In our example in Fig.|9|it is observed that we are not able to iden- 
tify foreground stars in a 10 mag Av cloud with a distance larger 
than 2 kpc when using colour excess extinction maps. Conversely, 
the star counting technique is able to detect this cloud at distances 
of up to 4 kpc at a 3 a level. For 3 mag Av the critical distance for 
the colour excess method is 3 kpc, compared to 4 kpc for the star 
counting technique. 

4.6 Suggested extinction determination method 

Given the results presented here, we suggest the following method 
to determine extinction maps: 

(i) Star count and colour excess extinction maps are calculated 
without foreground star corrections. These maps are used to deter- 
mine the opacity index /3, which is used to convert the colour excess 
maps into extinction maps. 

(ii) Foreground stars have to be selected. Extinction maps from 
colour excess are usually less noisy but are less reliable for distant 
clouds. Thus the extinction map that possesses the larger extinction 
values should be used for the foreground star selection. 

(iii) Now one needs to iterate the first two steps until no further 
foreground stars are found. Foreground stars cannot be identified in 
regions without extinction. We thus have to consider the accumu- 
lated star count diagram of the foreground stars when determining 
the extinction map using star counts. Using the colours of the iden- 
tified foreground stars and Monte-Carlo techniques, a colour offset 
for background regions can be determined. 

(iv) The final star count and colour excess maps are used to de- 
termine j3 and to convert the colour excess into extinction maps. 

(v) For a most accurate calibration determine the distance cor- 
rection using the cloud distance and coordinates, as well as the 
completeness limit of the observations. 



5 CONCLUSIONS 

• We have explained how a map of the near infrared extinction 
powerlaw index can be determined from colour excess, star count- 
ing, or a combination of both. A detailed analysis of the uncertain- 
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ties of each method is performed. We find that generally high S/N 
extinction measurements are required for a precision determination 
of (5. However, taking into account the cloud structure and possi- 
ble dependence of /3 on the column density, allows us to average 
different measurements and increase the accuracy. 

• We find that for 2MASS data of IC 1396 W the combination 
of star counting and colour excess methods lead to the highest sig- 
nal to noise /3-map. We find a dependence of the NIR extinction 
powerlaw index on the column density in the cloud. In low col- 
umn density regions /? is in agreement wi th the canonical value 
for the interstellar medium of 1.85 (Draine 1 2003)). At higher col- 
umn densities, however, /3 drops to values as low as 1.0. We fur- 
ther show that the column density profile of the cloud differs when 
measured by colour excess or star count extinction maps. Convert- 
ing the colour excess into extinction using the correct value of j3 
at each position, corrects these differences. It is thus vital to check 
and correct for any systematic changes of /? in dark clouds, prior 
to analysing the column density profile in extinction maps obtained 
from colour excess. 

• We present a scheme to test the accuracy of extinction deter- 
mination methods using artificial observational data. This is based 
on the model of stellar population synthesis of the Galaxy by Robin 
et al. «2003f) . We show that all extinction determination methods 
based on colour excess are subject to small but systematic offsets. 
These are due to the fact that the population of stars seen through a 
cloud is different from the population of stars in an extinction free 
control field (even if close by). The resulting offsets, however, can 
be determined with the knowledge of the cloud position, distance 
and extinction, as well as the completeness limit of the photomet- 
ric catalogue. The same scheme allows us to estimate the distance 
out to which clouds can be detected with star counts and colour ex- 
cess. Generally, star counts allow us to detect more distant clouds 
than colour excess. This is particularly useful for the first step of 
foreground star selection in the case of distant clouds. 
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